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ABSTRACT 



Aims. We present a catalog of point y-ray sources detected by the EGRET detector aboard the Compton Gamma Ray Observatory. 
We have used the whole y-ray dataset of reprocessed photons at energies above 100 MeV together with new Galactic interstellar 
emission models based on recent CO, HI, dark gas, and interstellar radiation field data. Two different assumptions have been used for 
the cosmic-ray distribution in the Galaxy to explore the resulting systematic uncertainties in source detection and characterization. 
Methods. We have used the same 2-dimensional maximum-likelihood detection method as for the 3rd EGRET catalogue. 
Results. The revised catalogue lists 188 sources, 14 of which are marked as confused, compared to the 271 entries of the 3rd EGRET 
(3EG) catalogue. 107 former sources have not been confirmed because of the additional structure in the interstellar background. The 
vast majority of them were unidentified and marked as possibly extended or confused in the 3EG catalogue. In particular, we do not 
confirm most of the 3EG sources associated with the local clouds of the Gould Belt. Alternatively, we find 30 new sources with no 
3EG counterpart. The new error circles for the confirmed 3EG sources largely overlap the previous ones, but several counterparts of 
particular interest that had been discussed in the litterature, such as Sgr A*, radiogalaxies and several microquasars are now found 
outside the error circles. We have cross-correlated the source positions with a large number of radio pulsars, pulsar wind nebulae, 
supernova remnants, OB associations, blazars and flat radiosources and we find a surprising large number of sources (87) at all 
latitudes with no counterpart among the potential y-ray emitters. 
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1. Introduction 

The Energetic Gamma-Ray Experiment Telescope (EGRET), 
which operated on board the Compton-Gamma Ray Observatory 
from April 1991 to May 2000, detected photons in the 20 MeV 
to 30 GeV range. The observation program made use of the 
large instrumental field of view (25° in radius) to cover the 
whole sky and for in-depth studies of specific regions. The re- 
sulting exposure and flux sensitivity to point sources are there- 
fore not uniform across the sky. The sensitivity threshold also 
varies because of the intense background emission that arises 
from cosmic-ray interactions with the interstellar gas and pho- 
ton fields in the Milky Way. The minimum flux that EGRET 
could detect steeply rises with decreasing Galactic latitude. In 
order to detect point sources and assess their significance in 
these varying conditions, a 2-dimensional maximum-likelihood 
method using binned maps had been developed for the COS- 
B d ata flPollock et al., 1981) and implemented for the EGRET 
one ( |Mattox et al., 1996| l. A first catalog using this method was 
published after 1.5 years of data (|Fichtel et al.,"l 994), followed 
by the second one (Thompson et al., 1995 1 and its supplement 



(Thompson et al., 1996 1 after 3 years of data. Lamb & Macomb 



(1997) presented a catalog of sources detected above 1 GeV. 
The last EGRET catalog (hereafter 3EG, |Hartman et al„ 1999| l 
comprised reprocessed data from April 1991 to October 1995 
with the interstellar emission model from Hunter et al. (1997) 
and extragalactic background from Sreekumar et al. (1998). 
This version contained 271 point sources including a solar flare, 
the Large Magellanic Cloud, five pulsars, one radiogalaxy de- 
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tection (Cen A), 66 high-confidence identifications of blazars 
(BL Lac objects and flat-spectrum radio quasars), and 27 lower- 
confidence blazar identifications. Because of the wide tails of 
the instrument point-spread function, seven potential artifacts 
were noted around the brightest sources and many sources were 
marked as confused or possibly extended. 

The 3EG catalogue also contained 170 sources with no at- 
tractive counterpart at lower energy. About 130 of them re- 
main unidentified as of today (see Grenier (2004) and refer- 
ences therein). Candidate counterparts that have been searched 
for include pulsars and their wind nebulae, supernova rem- 
nants, massive stars, X-ray binaries and microquasars, blazars 
and nearby radiogalaxies, luminous infrared and starburst galax- 
ies, and galaxy clusters. It was also noticed (Grenier, 1995 
|Grenier, 2 000, Geh rels et al., 2000] > that the most stable uniden- 
tified sources are significantly correlated with the nearby Gould 
Belt, a system of massive stars and interstellar clouds that sur- 
rounds the Sun at a distance of hundreds of parsecs. The offset 
position of the Sun with respect to the Belt centre and the Belt 
inclination of 17° to the Galactic plane indeed provides a useful 
spatial signature across the sky ( |Perrot & Grenier, 2003| ). 

EGRET went on observing for another 4.5 years after 
the 4 cycles used for the 3EG work. Its sensitivity was re- 
duced because of the ageing gas in the spark chamber, but 
it gathered nearly ten percent more photons and saw sev- 
eral new variable sources. Several authors ( |Nolan et al., 2 003 
|Sowards-Emmerd et al., 2005] ), however, noticed discrepancies 
between their studies and at least five 3EG sources. They 
failed to confirm sources and found others. The whole y- 
ray dataset and final instrument response functions have 
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also been significantly reprocessed by the EGRET team in 
2001. Furthermore, the spatial coverage of the CO surveys 
has reached higher latitudes since 1999, finding new small 
CO clouds ([Dame et al., 2001] >. In parallel, new HI surveys 
( |Kalberla et al., 2005) 1 have been completed to correct for the 
significant contamination of stray radiation in the older ones. 
Finally, an additional 'dark' gas component has been found in 
the Gould Belt clouds that significantly increases their mass and 
spatial extent (Gren ier et al., 2 005). The additional mass is struc- 
tured into large envelopes around the dense CO cores. They do 
not follow the HI and CO maps commonly used to trace atomic 
and molecular column-densities. So, the dark gas provides both 
y-ray intensity and structure that were not accounted for in the 
3EG background model. 

For all these reasons and in preparation of the new GLAST 
mission, it was necessary to revise the interstellar background 
model and to apply the EGRET detection method to the full nine 
years of data to build a new catalogue of sources above 100 MeV. 
In order to study the systematic uncertainties induced on source 
locations and fluxes by our limited knowledge of the intense in- 
terstellar background, we have applied the analysis to two dif- 
ferent background models exploiting the same new interstellar 
data, but using different approaches to constrain the cosmic -ray 
gradient across the Galaxy. 

2. The Galactic interstellar emission models 

The high-energy Galactic emission is produced by the interac- 
tion of energetic cosmic -ray electrons and protons with interstel- 
lar nucleons and photons. The decay of neutral pions produced 
in hadron collisions accounts for most of the emission above 300 
MeV. Inverse Compton (IC) scattering of the interstellar radia- 
tion field by electrons and their Bremsstrahlung emission in the 
interstellar gas are the other main contributors to the Galactic 
emission. The observed intensity therefore scales with the inte- 
gral along the line of sight of the cosmic-ray density times the 
gas or soft-photon one. 

The diffuse model used for the 3EG catalogue 
( |Hunter et al., 1 997 ) was based on a 3D-distribution of 
matter, cosmic -ray and soft-photon densities in the Galaxy, 
where the cosmic-ray density was assumed to be coupled to 
the gas one over a given length scale. This length as well as 
the CO-X.O-H2 conversion factor (X ratio) were adjusted to the 
data. The 3D gas map was obtained from the HI and CO line 
surveys and from kinematical distances derived for circular 
rotation. Distance ambiguities in the inner Galaxy were solved 
by splitting the gas into the far and near sides according to 
its expected scale height. Gas with velocities in excess of the 
tangent values was attributed to the tangent point and gas 
emission within 10° of the Galactic center and anticenter was 
interpolated from the regions just outside these boundaries and 
normalized to match the total emission seen along the line of 
sight. The resulting map is, however, still strongly biased to our 
side of the Galaxy, particularly for the atomic gas. This bias is 
reflected in the cosmic -ray density via the coupling length. 

For the present analyses, we have assumed an axisymmetric 
Galaxy for the cosmic -ray density and we have used gas column- 
density distributions in Galactocentric rings that are less subject 
to biases due to the strategy adopted to solve the cloud distance 
in the inner Galaxy. The radial velocity information in the HI 
and CO line surveys, together with the rotation curve of Clemens 
(1985) and the solar motion (v = 220 km/s at R = 8.5 kpc), have 
been used to partition the gas into 6 rings bounded by 3.5, 7.5, 
9.5, 1 1.5, and 13.5 kpc in Galactocentric distance (Digel et al., in 



preparation). Gas within 10° of the Galactic center and anticenter 
was interpolated as before. The all-sky Leiden-Argentina-Bonn 
(LAB) composite survey (Kalberla et al., 2005) was used for the 
HI data. Column densities, N(HI), were derived under the as- 
sumption of a constant spin temperature of 125 K. The velocity- 
integrated CO brightness temperature, W(CO), comes from the 
Center for Astrophysics compilation of observations at \b\ < 32° 
( |Dame et al., 2001) . The regions outside the survey boundaries 
should be free of bright CO emission. 

We have used two different approaches to account for the 
cosmic-ray density gradient. One is based on the Galprop model 
for cosmic-ray propagation developed by Strong et al. (2007, 
2004a, 2004b), using run number 49-6002029RB to derive the 
y-ray maps from pion decay, I,p, bremstrahlung radiation, hrem, 
and inverse Compton radiation, 7/c- This version includes sec- 
ondary electrons and positrons, an optimized cosmic-ray spec- 
trum to fit the GeV excess in the EGRET data, a cosmic -ray 
source distribution matching the radial profile of pulsars and su- 
pernova remnants, a radial gradient in the X factor, and the new 
HI and CO gas rings. 

The second model, hereby referred to as the Ring model, is 
based on the simpler, but realistic hypothesis that, if energetic 
cosmic rays uniformly penetrate all gas phases, the y-ray in- 
tensity in each direction can be modelled as a linear combina- 
tion of gas column-densities in the different rings, plus the IC 
intensity map (as predicted by Galprop), and an isotropic in- 
tensity (I/so) that accounts for very local IC emission and ex- 
tragalactic emission. This assumption has been used to derive 
gas emissivities in several rings from the COS-B and EGRET 
data ( |Strong et al., 1988| |Strong & Mattox, 1996) . We have re- 
produced these analyzes to derive gas emissivities for the new 
HI and CO rings using 9 years of EGRET data in three energy 
bands (> 100 MeV, 0.3 - 1 GeV, > 1 GeV). Both the Ring and 
Galprop models used the revised distribution of the interstellar 
radiation field ( |Porter & et al., 2005| |Moskalenko et al., 2006) 
to calculate the IC intensity map. The Galprop IC map is com- 
mon to both diffuse models. 

As indicated in the introduction, we have also included 
in the local ring the large column-densities of "dark" gas as- 
sociated with cold and anomalous dust at the transition be- 
tween the atomic and molecular phases (Gren ier et al., 2003) . 
This transitional phase is not traced in the radio. When remov- 
ing from total dust column-density maps the part that linearly 
correlates with N(HI) and W(CO), one is left with large en- 
velopes of excess dust around all the nearby CO clouds. The 
fact that the excess dust spatially correlates with significant dif- 
fuse gamma radiation indicates that cosmic rays pervade gas 
not accounted for in HI or CO. The gas-to-dust ratio in this 
phase, as inferred from the excess dust and correlated y-ray 
data, is normal. This phase appears to form an extended layer 
at the transition between the dense CO cores and the dens- 
est parts of the outer HI envelope of a cloud complex. It is 
best seen in total dust maps such as the reddening E(B-V) map 
( Schlegel et al., 1998), or low-frequency t hermal emission at 93 
GHz for WMAP (Finkbeine r et al., 1999) , or anomalous emis- 
sion near 20 GHz (Lagache, 2003). We constructed a "dark" gas 
column-density template, NHd ar k, by removing from the E(B- 
V) map the part linearly correlated with N(HI) and W(CO). 
This template was turned into gas column-densities by fitting 
it together with the N(HI) and W(CO) rings, as well as IC and 
isotropic components, to the all-sky y-ray maps. Because of its 
column-densities, dumpiness, and large spread across the sky 
(see Figure 4 in Grenier et al. (2005)), the "dark" gas component 
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Fig. 1. The top figure is the longitude profile of all photon counts observed by EGRET above 100 MeV at all latitudes (black error 
bars), compared with the diffuse counts predicted by the 3EG model (blue curve) and the Ring model (red curve). The bottom figure 
is the residual expressed in number of standard deviation, colors are the same as above, we added the Galprop residuals in purple. 
Counts from bright sources have been added to the diffuse component. For more visibility the plot are presented with a binning of 
4°. 



may strongly affect source detectability. This template was also 
added to the Galprop 49-6002029RB background model. 

To summarize, two diffuse backgrounds were constructed 
by fitting different components to the EGRET photon maps, in 
0.5° x 0.5° bins, in the three energy bands that will be used for 
source detection (> 100 MeV, 0.3 - 1 GeV, > 1 GeV). 

1 . With the Ring model, the predicted count rates are calculated 

as: 

Npredd, b) = [brings qHI,iN H l(n, I, b) + brings qCO,iW C o(n, I, 

+qdarkNH dark (l, b) + qichciU b) + I iso ] x e(l, b) 

+ 2,j=source S b j) fj PS F ( l h b j) 

2. and the Galprop model as: 

N pred (l, b) = [qjfiljfiil b) + qbremhremd b) + qdarkNH dark (l, b) 

+qichc(U b) + I Lm ] x e(l, b) 

e(lj,bj)fjPSF(lj,bj) 

In both models, e(l, b) and fj note the EGRET exposure map 
and source fluxes. The diffuse maps times the exposure were 
convolved with the EGRET PSF for an input E~ 2A spectrum be- 
fore adding the source maps. The EGRET count and exposure 
maps, the 3EG diffuse model, as well as the latest instrument 
response functions, were downloaded from the CGRO Science 
Support Center. They differ from those used for 3EG since 
they were reprocessed in 2001. The q parameters (gas emissiv- 
ities or relative contributions of different radiation components) 



were fitted to the data by means of a maximum likelihood with 
Poisson statistics. To avoid biasing the interstellar parameters, 
the model included the brightest sources detected during a first 
source detection iteration with a significance > 5cr, with fixed 
fluxes. Changing these fluxes within their statistical uncertain- 
ties do not significantly change the diffuse results. 

The resulting emissivities corresponding to the local gas are 
-fully consistant with Grenier et al. (2005) Table 1 . The emissivity 
gradient in the Galactic plane will be described in a separate pa- 
per. The quality of the fit can be seen in Figure Q] The top figure 
displays the longitude profile of all the EGRET photon counts 
above 100 MeV. The error bars are only statistical. The plot com- 
pares the best fit that can be obtained using the former 3EG dif- 

tise model with the longitude profile resulting from the present 
ing model. The bottom plot shows the longitude profile of the 
residuals and the improvement of the ring model over the 3EG 
one. It also shows the residuals for the best fit Galprop model. 
All modelled profiles include the brightest sources. Systematic 
differences can be seen in various places where the 3EG model 
significantly over-predicts and under-predicts the data while the 
new models behave better. Because of its larger flexibility (the 
gas emissivity gradient due to cosmic -ray variations is measured, 
not inferred from propagation properties or gas coupling), the 
Ring model was found to best fit the data. It is worth noting than 
even if the agreement is excellent, there still exists small devia- 
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Fig. 2. Map in Galactic coordinates of the residuals (expressed in 
o" = y/Npred values) between the E > 100 MeV photon counts 
(in 0.5° bin) and the best fit with the Ring model using Equation 



tions that can significantly impact source detection and charac- 
terization. 

The residual count map obtained above 100 MeV with the 
Ring model is presented in Figure 2. It displays the statistical dif- 
ference (N b s - N pre d)l JNp re d between the observed counts and 
those predicted from the diffuse background and bright sources 
using equation[T] The model globally fits very well the data. The 
extended blue fan-like structures with negative residuals are cor- 
related with the edge of several observing periods. They prob- 
ably result from a wrong exposure estimate at large angle from 
the instrument axis. They are visible independently of the choice 
of diffuse model (Ring, Galprop, or 3EG). Their spatial extent is 
large enough compared to the PSF size not to severely affect 
source detection, yet source fluxes in these directions are under- 
estimated. Uncertain knowledge of the off-axis instrument ex- 
posure is also reflected in the small model deficit (orange edge) 
bordering the fan-like excesses. We have checked for suspicious 
strings of faint sources that would correlate with these instru- 
mental features. 

The use of two different background models allowed us 
to study their impact on source detection and characterization. 
Given its higher likelihood value and locally flatter residuals, 
the Ring model was used to derive the default source flux and 
location. The values obtained with the Galprop background are 
used to illustrate the amplitude of the systematic uncertainty due 
to the background modelling. When searching for sources we 
used the diffuse emission parameters calculated from this global 
fit. We adjusted a source flux together with a free normaliza- 
tion of the total diffuse flux within 15° around each pixel, and a 
free isotropic flux. This procedure is the same as used for 3EG 
(Gmult and Gbias). These two parameters correct for small lo- 
cal mismatches between the diffuse model and the data. Gmult 
fluctuates around 1 . 



3. Source detection 

As for the derivation of the 3EG catalogue, we have used the 
LIKE code (Mattox, 1996, version 5.61) to compute the 2- 
dimension binned Poisson likelihood of detecting a source at a 
particular location on top of the diffuse background. LIKE cal- 



culates the Test Statistic (TS) value that compares the likelihood 
of detecting a PSF-like excess above the background to the null 
hypothesis - a random background fluctuation - for a given po- 
sition. The likelihood (L,) is calculated as the product, for all 
pixels within 15° of a specific position, of the Poisson probabil- 
ities of observing photons in a pixel where the number of counts 
is predicted by the model (background + source). The likelihood 
ratio test statistic is defined as TS = -KLtiLq — LnL\), where 
the likelihood values L\ and Lq are respectively optimized with 
and without a source in the model. Asymptotically, the TS dis- 
tribution follows ax 2 one. The detection significance of a source 
at the given position is y/TScr (Mattox 1996). 

Sources have been searched for in the summed maps corre- 
sponding to cycle 1, 2, 3, 4, 1+2, 3+4, 1+2+3+4, 5+6, 7+8+9, 
1+2+3+4+5+6+7+8+9. In addition, we have analyzed the 46 
individual periods listed in Table[TJfor which flaring 3EG sources 
had been detected. As for the summed maps, the individual pe- 
riod maps retained only photons with inclinations within 30° 
from the instrument axis, or 19° for cycle 6, 7, 8, and 9. Photons 
and exposure maps were binned to 0.5° x 0.5°. 

To build the 3EG catalogue, sources were detected only in 
the integrated E > 100 MeV band. TS maps were then con- 
structed in three energy bands (> 100 MeV, 0.3 - 1 GeV, and > 1 
GeV) from the observation (single or summed) with highest TS 
and a source final position was obtained from the smallest er- 
ror contours. Given the modern computer performance, we have 
directly searched for sources independently in the three energy 
bands. 

At 100 MeV, the EGRET PSF is wide and there exists dis- 
crepancies between its real shape, as observed in bright sources, 
and the modelled one. In practice, differences may also come 
from a more complex source spectrum than the single power- 
law assumed to integrate the PSF. A choice of 300 MeV instead 
of 100 MeV for the lower analysis threshold might have been a 
better trade-of between count rates for detection and systematic 
uncertainties in the PSF. We have, however, kept a lower limit 
of 100 MeV as in 3EG in order to account for soft sources and 
to allow comparison with the 3EG results. We have assumed a 
spectral index of 2.0 for all sources but for 1 1 bright ones which 
had a 3EG spectral index far from 2.0. For the latter, we have 
used their 3EG index to integrate the PSF. 

Each of the 10 all-sky summed maps was divided, both in 
Galactic and equatorial coordinates, in 45 zones with a large 
overlap. The use of both coordinates systems is required since 
source images are deformed in rectangular projection at high 
latitude or declination. For each zone, each individual period, 
and each of the 3 energy bands (> 100 MeV, 0.3 - 1 GeV, and 
> 1 GeV), we calculated a TS map for excesses above the back- 
ground. Sources were iteratively detected from high TS to low 
TS in successive TS maps. Between each steps, the detected 
sources were included in the background model until no ex- 
cess with y/TS > 3 was left in the final TS map. An exam- 
ple of the iteration around Geminga is given in Figure 3. Peaks 
in the TS map were automatically detected with SExtractor 
( |Bertin & Arno uts, 1996) and converted into source position by 
taking the TS -weighted centroid in the region enclosed by the 
95% confidence contour around this position. Source positions 
were recalculated at each iteration to take into account the in- 
fluence of the neighbouring sources. More than 1100 rS-maps 
were thus calculated at the CCIN2P3 Computing Center. 
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Name 

2+ 

0020 

0040 

0050 

0200 

0210 

0220 

0230 

0250 

0260 

0290 

36+ 

0420 

0430 

0440 



Sum of viewing periods 
0002+0003+0004+0005 



0360+0365 



Name 

2040 

virg2 

2110 

2230 

2260 

227+ 

229+ 

2310 

3023 

314+ 

3170 

319+ 

3200 

328+ 

3290 



Sum of viewing periods 



2040+2050+2060 



2270+2280 
2290+2295 



3140+3150 



3190+3195 



3280+3310+3315+3330 



Name 

3315 

330+ 

335+ 

vrg3a 

3355 

3360 

3385 

3390 

4040 

4100 

4130 

4180 

419+ 

4210 

4230 

4235 



Sum of viewing periods 

3300+3320 
3350+3355 

3040+3050+3060+3070+3080+3086 



4191+4195 
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Fig. 3. An example of the iterative source detection with the 2D 
binned likelihood around Geminga at energies above 100 MeV. 4 
consecutive TS maps are shown. Sources are detected, then are 
included in the background for the next step until no significant 
one is left. The colourbar gives TS . 

4. Catalogue construction 

To account for real versus modelled PSF discrepancies in ex- 
tremely bright sources, for instance to account for the splitting 
in two of the bright pulsar sources or for the artifacts in the 
Vela tails, we have removed all the source candidates within 3.5° 
of the intense sources (that exhibit more than 800 photons in a 
map). For less intense sources, we have checked the probability 
of having a double versus single source with a specific likeli- 



hood calculation, using the likelihood ratio between the 2 cases 
to keep or reject the double source. 

At the end of this stage, most sources have two possible po- 
sitions per energy band and observation, one from the Galactic 
coordinate map and one from the equatorial one. We cross- 
compared the two and selected the position from the least de- 
formed projection. Sources detected only once were not included 
in the list unless their latitude or declination were higher than 40° 
or their longitude or right-ascension were less than 5° from the 
map edges. 

At this stage, most sources have three possible positions 
(with energy) for a given observation. We chose among the three 
the position corresponding to the smallest 95% confidence con- 
tour, unless its peak ^TS were 1.5 smaller than found in an 
another energy band. The latter condition reduces the risk of in- 
correct source assignment during the cross-comparison phase. 
Sources found at low energy, but not at high energy were in- 
cluded in the list, as well as sources found only at high energy. 

We have used the same criteria to cross-compare the source 
positions for individual periods and summed cycles in order to 
obtain a final list of candidate sources with the best position 
from the different energy bands and periods/cycles. We followed 
the whole procedure with both the Ring and Galprop interstellar 
backgrounds. We obtained respectively 1192 and 1225 candi- 
date sources with the Ring and Galprop models. Source fluxes 
and y/TS values above 100 MeV were calculated for these sets 
of positions for the different periods and cycles. Unlike in 3EG, 
we did not adjust the position of the identified sources (AGN or 
pulsars) to that of their radio counterpart. 

We adopted the same detection threshold as for the 3EG cat- 
alogue ( y/TS > 5 at \b\ < 10° and y/fS > 4 elsewhere) and 
found 188 and 208 significant sources for the Ring and Galprop 
models, respectively. We manually checked the TS maps of all 
the sources that barely passed the detection threshold with the 
Ring model and had y/TS ~ 3 with the Galprop one. 

We emphasize the fact that the order and criteria applied to 
cross-correlate positions between the excesses detected in differ- 
ent energy bands and time periods can strongly affect the catalog 
list near the detection threshold. Several strategies were tested 
before adopting the present one, but one must remember that 
a faint source can pass or drop below the threshold by slightly 
changing its position or that of its neighbours. Given the steep 
increase in source numbers with decreasing TS , we also empha- 
size that a small change in the TS threshold, alternatively in the 
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background over which the source TS is calculated, results in 
a large change in the number of catalogue entries. For instance, 
lowering the y/TS threshold by 0. 1 would add 27 sources. 

5. Catalogue description 

The EGR acronym has been adopted for the EGret Revised 
source list presented in Tableland Table lA. fl in a format similar 
to the 3EG one. As explained above, the source characteristics 
(position and flux, and their uncertainties) have been determined 
with the Ring model because of its higher flexibility, better fit, 
and flatter residual map. A secondary position and flux has been 
measured with the Galprop model and is listed in Table [2] and 
Table lA. f| to illustrate the amplitude of the systematic uncertain- 
ties due to the choice of interstellar model. 

Sources found within a radius of 1 .5 PSF FWHM from a very 
bright source, and/or with very asymmetric TS map contours are 
not included in Tableland Table |A~T1 Still, they represent sig- 
nificant excesses of photons above the background which may 
be due to extended sources, or structures not properly modelled 
in the interstellar emission, or artifacts due to incorrect PSF tails. 
This list of 14 confused sources is given in Table IB. fl under the 
acronym EGRc for EGret Revised confused. 

For both tables, the description for each column follows: 

1 . Num: source number in order of increasing right ascension. 

2. Name: source name based on J2000 coordinates. 

3. RA and Dec: J2000 equatorial coordinates in degrees. 

4. 1 and b: Galactic coordinates in degrees. 

5. 695 : angular radius, in degrees, of a circular cone which con- 
tains the same solid angle as the 95% confidence contour. 

6. F: flux in 1(T 8 photon cirT 2 s -1 for E > 100 MeV and for 
each time period. 

7. cr F : lcr statistical flux uncertainty in 10~ 8 photon cirT 2 s _1 . 



9. 
10. 



11. 



12. 



13. 



Cnts: number of photons detected with E > 100 MeV. 

y/TS : statistical significance of the detection. 

vp: short viewing period as defined in Table Q] or summed 

cycles noted px for cycle x, pijkl for the sum of cycles i, j, 

k, and 1, and p\9 for the total of 9 cycles. 

\ sys and b JVJ : Galactic longitude and latitude obtained with 

the Galprop background model. 

F !V! : flux obtained with the Galprop background model, in 
10 -8 photon cirT 2 s . 

3EG: third EGRET catalog counterpart source name if one 
exists within a radius of 1 PSF FWHM (2° for E > 100 MeV) 
from the EGR source and if the nearest neighbour relation 
between the EGR and 3EG sources is univocal (the nearest 
neighbour of the EGR source is the 3EG one and vice versa). 
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Fig. 4. Spatial distribution, in Galactic coordinates, of the EGR 
sources. The confused sources are marked as open circles. 
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Fig. 5. Spatial distribution, in Galactic coordinates, of the 3EG 
sources. 
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Fig. 6. Spatial distribution, in Galactic coordinates, of the 3EG 
sources with no counterpart in EGR: the unidentified sources as 
circles and the identified AGN as stars. The filled circles and 
stars mark the sources that were flagged as extended or confused 
in the 3EG catalogue. 



6. Comparison with the 3EG catalogue 

The revised catalogue contains 174 sources plus 14 confused 
sources compared to the 265 entries of the 3EG catalogue (ex- 
cluding the Vela artifacts). Their spatial distribution across the 
sky looks different from that of the 3EG sources, as illustrated in 
Figures |4] and |5] The accumulation of faint 3EG sources within 
30° of the Galactic center is much more reduced in the new re- 
sults and fewer sources are seen below 30° in general. These 
changes at low and mid latitudes are primarily due to the in- 
crease in background intensity from new HI, CO, and dark gas 
structures. At high latitude, the use of more y-ray observations 
and of a revised large-scale IC component in the background 
may also explain why a handful of 3EG sources have fallen be- 
low the detection threshold whereas new ones are now detected. 



The names of the 107 unconfirmed 3EG sources are listed in 
Tableland they are displayed in Figure [6] They comprise only 
six sources that had been firmly identified as AGN by Hartman 
et al. (1999), but that had been flagged as extended or con- 
fused by the EGRET team. In fact, the proportion of these ex- 
tended or confused cases among the unconfirmed 3EG sources 
is overwhelming (95%) and significantly larger than among the 
confirmed ones. The unconfirmed and confirmed 3EG groups 
respectively show 69% and 33% of possibly extended 'em' 
sources. Figure [6] also shows that the vast majority of uncon- 
firmed 3EG sources were unidentified and spatially correlated 
with the Gould Belt system of nearby clouds. They follow the 
characteristic trace of the inclined Belt across the sky, gathering 
at \b\ < 30°, more at positive latitudes toward the Galactic centre 
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Num Name 


RA 


Dec 


1 


b 


095 


F 


0"F 


Cnts 




vp 


1st, 




F„., 


3EG 


1 EGR J0O08+73OS 


2.01 


73.14 


119.75 


10.54 


0.20 


39.7 


4.4 


330 


10.9 


pl9 


119.75 


10.54 


41.0 


3EGJ0010+7309 














63.9 


11.6 


96 


7.2 


Pi 






















33.4 


9.6 


61 


4.1 


P2 






















22.4 


8.7 


37 


3.0 


P4 






















48.8 


7.4 


162 


8.2 


p 12 






















21.6 


7.3 


52 


3.4 


p34 






















37.0 


5.3 


212 


8.5 


pi 234 






















44.6 


8.3 


1 15 


6.6 


p56 






















33.1 


9.7 


60 


4.0 


2110 










2 EGR J0028+0457 


7.06 


4.95 


112.15 


-57.44 


0.51 


14.3 


4.6 


31 


4.1 


p34 


112.15 


-57.44 


14.4 
















7.2 


4.8 


13 


1.7 


Pi 






















13.9 


5.9 


20 


3.0 


p3 






















14.0 


7.2 


10 


2.7 


p4 






















7.2 


4.8 


13 


1.7 


pl2 






















10.7 


3.3 


43 


4.0 


pi 234 






















10.4 


3. 1 


46 


4. 1 


pl9 






















24.2 


11.1 


14 


2.9 


3200 






















24.7 


15.6 


6 


2.3 


3360 










3 EGR J0039-0945 


9.75 


-9.75 


112.76 


-72.38 


0.27 


13.0 


3.5 


48 


4.8 


pl9 


112.65 


-72.40 


13.1 


3EGJ0038-0949 














14.6 


5.8 


23 


3.4 


Pi 






















15.7 


5.4 


24 


4.0 


p4 






















14.6 


5.8 


23 


3.4 


pl2 






















11.0 


4.5 


22 


3.2 


p34 






















12.3 


3.6 


43 


4.5 


pi 234 






















22.2 


17.7 


5 


1.6 


p789 











new EGR sources 



Table 3. Names of the 3EG sources with no EGR counterpart 




Fig. 7. Spatial distribution, in Galactic coordinates, of the new 
EGR sources with no 3EG counterpart. The confused sources 
are marked as open circles. 



and below the plane in the anticenter. The EGR source sky dis- 
tribution in Figure [4] does not exhibit the Gould Belt signature 
anymore. 

The fact that many 3EG sources are not confirmed in the 
present analyses should not cast doubts on the detection method 
from a statistical point of view. They did correspond to signif- 
icant photon excesses above the background in the 3EG analy- 
ses, but, in the absence of some structures in the predicted inter- 
stellar background, an ensemble of point sources with the wide 
EGRET PSF would compensate for the missing clouds and yield 
an excellent fit to the data. Figure[8]illustrates this point with the 
unidentified source 3EG J0556+0409 detected at 7.2 cr in 3EG. 
The left side shows the TS-map corresponding to the second 
stage of the iterative source detection around Geminga above 
100 MeV. It is the same as in Figure[3]but we have used here the 
3EG diffuse emission model instead of the Ring one. The same 
sources are detected except for 3EG J0556+0409 which is not 
seen in Figure [3] Instead an excess of diffuse emission appears 
in the ratio of the Ring to 3EG background intensities (Figure 
[8j right). The photons attributed to a point source in 3EG where 
in fact coming from a gas cloud in the Galaxy. This is prob- 
ably still the case in the present analysis, although to a lesser 
degree, in particular at very low latitude where optical thickness 



3EG JO 130- 175 8 
3EG J0348+3510 
3EG J0416+3650 
3EG J0439+1555 
3EG J0459+0544 
3EG J05 10+5545 
3EG J0533+4751 
3EG J0546+3948 
3EG J0622-1139 
3EG J0702-6212 
3EG J0808-5344 
3EG J1013-5915 
3EG J1052+5718 
3EG J 1227+4302 
3EG J 1300-4406 
3EG J1316-5244 
3EG J 1329-4602 
3EG J1527-2358 
3EG J 1627-24 19 
3EG J1633-3216 
3EG J 1639-4702 
3EG J1653-2133 
3EG J1709-0828 
3EG J1718-3313 
3EG J1735-1500 
3EG J 1744-03 10 
3EG J1757-0711 
3EG J1810-1032 
3EG J1824-1514 
3EG J 1834-2803 
3EG J 1850-2652 
3EG J 1904- 1124 
3EG J 195 8-4443 
3EG J2022+4317 
3EG J2 100+60 12 
3EG J2255+1943 



3EG J0245+1758 
3EG J0404+0700 
3EG J0426+1333 
3EG J0439+1105 
3EG J0459+3352 
3EG J0520+2556 
3EG J0542+2610 
3EG J0556+0409 
3EG J0628+1847 
3EG J0706-3837 
3EG J0821-5814 
3EG J1014-5705 
3EG J1212+2304 
3EG J1235+0233 
3EG J 1308+8744 
3EG J 1323+2200 
3EG J 1447-3936 
3EG J1600-0351 
3EG J1631-1018 
3EG J 1634- 1434 
3EG J 1646-0704 
3EG J1659-6251 
3EG J1714-3857 
3EG J1720-7820 
3EG J 174 1-2050 
3EG J 1744-30 11 
3EG J1800-0146 
3EG J1823-1314 
3EG J1825+2854 
3EG J1836-4933 
3EG J1858-2137 
3EG J1928+1733 
3EG J2016+3657 
3EG J2034-3110 
3EG J2206+6602 
3EG J2359+2041 



3EG J0323+5122 
3EG J0407+1710 
3EG J0435+6137 
3EG J0458-4635 
3EG J0500+2529 
3EG J0521+2147 
3EG J0542-0655 
3EG J06 16-0720 
3EG J0634+0521 
3EG J0747-3412 
3EG J09 10+6556 
3EG J1045-7630 
3EG J1222+2315 
3EG J 1249-8330 
3EG J1308-6112 
3EG J1329+1708 
3EG J 1500-3509 
3EG J 16 16-2221 
3EGJ163 1-4033 
3EG J1635-1751 
3EG J1649-1611 
3EG J 1704-4732 
3EGJ17 17-2737 
3EG J1733+6017 
3EG J1741-2312 
3EG J1744-3934 
3EG J1806-5005 
3EG J 1824+ 3441 
3EG J1828+0142 
3EG J 1850+5903 
3EG J 1903+0550 
3EG J1958+2909 
3EG J2020-1545 
3EG J2035+4441 
3EG J22 19-7941 



in HI and CO severely limits our knowledge of the true column- 
densities. Other sources may also be due to increased cosmic -ray 
densities in specific clouds with respect to the local Galactic av- 
erage. Over-irradiated clouds near cosmic-ray sources would be 
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Table 4. Names of the new EGR sources with no 3EG counter- 
part 



EGR J0028+0457 
EGR J0141 + 1719 
EGR J0509+0550 
EGR Jl 158-1950 
EGR J1642+3940 
EGR J 1920+4625 
EGR J2202+3340 
EGR J2308+3645 
EGRc J09 12+7 146 
EGRc J 1255-0404 



EGR J0057-7839 
EGR J0243-5930 
EGR J0540+0657 
EGRJ1259-2209 
EGR J 1740+4946 
EGRJ1959+4322 
EGR J2233-4812 
EGRc J08 18-46 13 
EGRc J0927+6054 
EGRc J1332-1217 



EGR JO 100+4927 
EGR J04 13-3742 
EGR J 1122-5946 
EGR J1619+2223 
EGR J1814+2932 
EGR J2027-4206 
EGR J225 8-2745 
EGRc J0842-4501 
EGRc J1038-5724 
EGRc J22 15 +065 3 




Fig. 8. Second stage of the iterative source detection around 
Geminga (see Figure obtained using the 3EG model (left) 
and map of the Ring model intensity divided by the 3EG one 
(right). The excess in the TS map assigned in 3EG to the 3EG 
J0556+0409 point source corresponds to a local underestimation 
of the diffuse emission in the 3EG model. Maps are given in 0.5° 
bins and galactic coordinates 



3EG 95 



Fig. 10. Histogram of the relative angular separation between the 
positions found for the EGR and 3EG counterparts in units of the 
95% confidence angle for each source. 




l F EOR- F s y s,l 7CT F 

Fig. 11. Histogram of the relative flux differences \Fegr ~ 
F sys \l<TF EGR measured with the Ring and Galprop models in units 
of the statistical error on flux for each source. All fluxes are mea- 
sured above 100 MeV. 




Fig. 9. Histogram of the relative flux differences \Fecr ~ 
FiEc\/o'F EGR measured between the EGR and 3EG counterparts 
in units of the statistical error on flux for each source. All fluxes 
are measured above 100 MeV. 



detected as a single or cluster of point sources, depending on 
their angular scale. 

For the 8 1 EGR sources that do have a 3EG counterpart, we 
find a reasonable agreement in position and flux from both anal- 
yses. On average, we find 3% lower fluxes in the EGR analysis 
with respect to the 3EG one because of the increase in Galactic 
background. Figure [9] shows the histogram of ratios of the EGR 
and 3EG flux difference over the statistical error on flux for each 



source: \Fecr - Fj,ec\Io'f egr - The EGR flux was taken for the 
observation with highest y/TS and compared to the 3EG coun- 
terpart flux for the same time period if available. Average P19 
fluxes were compared to the 3EG PI 234 average for non flaring 
sources. The flux differences are modest (17% rms dispersion) 
and in most cases smaller than the statistical uncertainties on 
flux estimates. Similarly, FigurefTUlshows that the angular sepa- 
rations between EGR and 3EG counterparts are often consistent 
with the 095 error radii. Yet, thirty sources have been found as 
far as 0.5° from the 3EG position and this will greatly impact 
counterpart searches and identification at other wavelengths. 

On the other hand, we find 30 new EGR sources with no 
3EG counterpart. Their names are listed in Tableland they are 
displayed in Figure [7] Most of them are detected just above the 
threshold and 1 1 of them were indeed present in the 3EG com- 
plementary list, just below the significance threshold. 

7. EGR source distributions and potential 
counterparts 

Because of the new gas data we have used at intermediate lati- 
tude, the comparison between the EGR and 3EG source charac- 
teristics allows to judge, to some extent, the impact of our limited 
knowledge of gas mass tracers. The comparison between the flux 
and positions obtained with the Ring and Galprop models gives 
an estimate of the systematic uncertainties due to our limited 
knowledge of the true cosmic -ray distribution across the Galaxy. 
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Fig. 12. Histogram of the relative angular separation between the 
positions found with the Ring and Galprop models in units of the 
95% confidence angle for each source. 




40 60 

latitude (deg) 

Fig. 14. Latitude distribution of the EGR sources with young 
Galactic sources at \b\ < 3°, nearly isotropically distributed 
sources far from the plane, as expected from the black curve, 
and a flattening at mid-latitude because of the rapid increase in 
the interstellar background flux. 



FigureQ~T]and FigureQ~2]show that, in most cases, the differences 
are smaller than the statistical uncertainties. The distribution of 
95% confidence radii peaks between ~ 0.2° and ~ 0.7°. The 
uncertainty in the background induces an additional systematic 
error of ~ 0.2° for most sources. It should be kept in mind while 
looking for counterparts. 

We have searched the EGR error circles for poten- 
tial counterparts of interest such as pulsars from the 
ATNF catalogue (Man chester et al., 2005|> , blazar candidates 
from the ASDC list (|Massaro et al., 2005| l and the CGRaBS 
list ( |Healey et al., 2008| l, other flat radio sources from the 
CRATES compilation (Healey et al., 2007 1, supernova rem- 
nants from the Green catalogue (Green T2006| l, OB associations 
( |MerNik&Efremov(1995)|l, and X-ray an d TeV pulsar wind 
nebulae (Li et al. (2008) and Grenier, 20 08). The results are dis- 
played in Figure 1131 We have found 13 radio pulsar associa- 
tions in addition to the 6 objects firmly identified by EGRET. 
Thirteen EGR sources coincide with supernova remnants, 9 with 
pulsar wind nebulae, 7 with OB associations, 53 with blazar 
candidates, and 19 with other flat radiosources. These associa- 
tions should not be considered as identification, but as spatial 
coincidences worthy of further investigations, in particular with 
the improved angular resolution of GLAST. Yet, they reveal that 
as many as 87 sources have no obvious counterpart among the 
well-known y-ray emitters despite the large number of pulsars 
(1775) and radiosources (11 000) that were cross-correlated with 
the sources and that spread across the entire sky and along the 
Galactic plane. The lack of blazar counterparts is all the more 
surprising that the spatial distribution of the sources off the plane 
is quite reminiscent of an isotropic, therefore extragalactic, dis- 
tribution. The latitude distribution, shown in Figure [14] is quite 
consistent above 30° with a sample drawn from a uniform pop- 
ulation according to the exposure map, as shown by the black 
curve. The distribution flattens at lower latitude because of the 
increased background that drastically limits the survey sensitiv- 
ity. Studying the consistency with an extragalactic population 
at medium latitudes and the implication of the lack of a flat ra- 
diosource is beyond the scope of this paper and will be addressed 
in a separate one. The sharp peak below 3° in latitude indicates 
young emitters. Their clustering in the inner Galaxy (I < 30°), 
toward the direction tangent to the Carina arm, and toward the 
Cygnus region outlines their close relationship to large molec- 
ular complexes and star forming regions at a distance of a few 
kpc. 




Fig. 15. TS -map obtained at energies above 1 GeV toward the 
Galactic center above the 3EG (a) and Ring (b) interstellar model 



8. Discussion on specific sources 

There is considerable interest in the physical processes occurring 
in the Galactic center region. The 3EG catalogue lists one source 
located toward the Galactic center, 3EG J1744-3011. We find 
two point sources in this region, EGR J 1740-2851 at I = -0.55°, 
b = 1.05° and EGR J1747-2852 at I = 0.21°, b - -0.24°. Figure 
[15] display the rS-map for photons with energies above 1 GeV 
above the 3EG and the Ring background models. The #95 error 
radius around EGR J 1740-2851 and EGR J 1747-2852 formerly 
excludes the Galactic Center but source locations and fluxes in 
this direction should be taken with extreme caution since the 
high gas optical depth around the Galactic center and the ve- 
locity pile-up toward the center induce large uncertainties in the 
total gas column densities. 

Coincidences with supernova remnants were noted 
(Sturner & Dermer(1995)) and are confirmed in the present 
analysis (see Table |5J, but several also host a pulsar wind 
nebula, as in CTA 1 and IC 443, so we need much higher 
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Fig. 13. The revised EGRET source catalog, shown in Galactic coordinates. The symbols indicate the counterpart types found in 
the error box: identified pulsars as black squares; other ATNF pulsars as open squares; LSI +61 303, LMC, and solar flare as 
black triangles; ASDC and CGRaBS blazar candidates as black diamonds; other flat-spectrum radio-sources from CRATES as open 
diamonds; supernova remnants from the Green catalogue as stars; no counterpart as crosses. 



Table 5. Names of the sources and supernova remnants found in 
spatial coincidence 



EGRJ0008+7308 


G119.5+10.2 


CTA1 


EGRJ0617+2238 


G189.1+3.0 


IC443 


EGRJ0633+0646 


G205.5+0.5 


Monoceros 


EGRJ1710-4435 


G343.0-6.0 


RCW114 


EGRJ 1800-2328 


G6.4-0.1 


W28 


EGRJ 1800-2328 


G6.5-0.4 




EGRJ1838-0420 


G27. 8+0.6 




EGRJ1838-0420 


G28.8+1.5 




EGRJ2020+4019 


G78.2+2.1 


y Cygni 


EGRJ2227+6114 


G106.3+2.7 




EGRJ0225+6240 


G132.7+1.3 


HB3 



PSR J0030+0451 



riqht ascension 



resolution y-ray images to identify the origin of the emission, 
especially in these crowded regions. EGRET detections are 
confirmed toward two TeV-emitting wind nebulae around 
PSR J1420-6048 (in Kookaburra, EGRJ 14 18-6040) and PSR 
J1826-1334 (EGRJ1825-1325). Another interesting candidate is 
the wind nebula of the 1 1-kyr old and very energetic pulsar PSR 
J2229+61 14 toward EGRJ2227+61 14. 

We also note, as shown in Figure [16] the positional coinci- 
dence within 0.5° between the new EGR J0028+0457 source and 
the millisecond X-ray pulsar PS R J0030+045 1 This 300 pc dis - 
tant pulsar, discovered in 2000 ( |Somer, 2000) ID'Amico, 20001 , 
has an X-ray counterpart exhibiting a double peaked pulse pro- 
file as seen by ROSAT (Beck er et al., 20 00). Millisecond pulsars 
have low magnetic fields, they produce relatively few electron- 
positron pairs so the electric field is not screened and the spectral 
cutoff due to pair production attenuation occurs at high energy. 
They are therefore good candidate for accelerating particles to 
high energies. Harding et al. (2005) has predicted a y-ray flux 
for PSR J0030+045 1 well above the one of the y-ray millisec- 



Fig. 16. Likelihood TS contours for energies above 100 MeV 
and periods incompassing PSR J0030+045 1 . The cross, the plus 
sign and the black dot respectively mark the EGR catalog posi- 
tion, the position with maximum likelihood and the pulsar loca- 
tion 



ond pulsar PSR J0218+4232 for which a pulsed emission was 
marginally detected (Kuiper et al., 2000 1. 

Four massive binaries have been detected at TeV energies, 
namely PSR B 1259 -63 dAharonian et al., 2005), LSI +61° 303 
dAlbert et al., 2006)1, LS 5039 ( |Aharonian et al., 2006)) , and Cyg 
X-l ( ] Albert et al., 20 07), thus illustrating very efficient particle 
acceleration in compressed or shocked pulsar winds, as well as 
in microquasar jets. Inverse Compton scattering of the bright 
stellar radiation would dominate at GeV energies. We find no 
interesting EGRET counterpart to these high-energy objects, but 
for the LSI +61° 303 radiosource. The latter had long been as- 
sociated with the COS-B source 2CG 135+01 and the EGRET 
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65 64 63 62 61 

galactic longitude 



Fig. 17. Likelihood TS contours ( 50%, 68%, 95% and 99% con- 
fidence) for energies above 100 MeV and period 5 190. The cross 
is the EGR catalog position, the plus sign the position with max- 
imum likelihood and the black dots mark the radio positions of 
3C345, Mrk 501, and 4C+38.41 



source 2EG J0241+61 19 ( |Kniffen et al., 1997| l, yet it had moved 
out of the 3EG error box and the marginal y-ray variability could 
not be associated with the radio flux variations. In the present 
analysis, we find the radiosource very near the centre of the EGR 
J0240+6112 source. On the other hand, we find no source to- 
ward the dust enshrouded microquasar candidate, AX J1639.0- 
4642, or the Be/X-ray binary, AO 0535+26, both proposed as 
3EG counterparts ( |Combi et al., 2003[|Romero et al., 200T] |. 

Another noticeable new source is EGR J 1642+3940 de- 
tected at 5.8cr rather close to 3C345. 3C345 is one of the most 
prominent flat spectrum (a = -0.1) radio-loud, superluminal 
sources and is therefore an excellent candidate for a y-ray blazar. 
EGRET has viewed this region 12 times, in particular during 
period 5190 when a flare was found. We have analyzed again 
this particular period with the Ring model since it had not been 
used in the overall detection search. Figure [17] shows the result- 
ing TS contour for photons with energies above 100 MeV that 
is well centered on 3C345. The cross corresponds to the EGR 
position (period 5190), the plus sign to the position with maxi- 
mum likelihood and the black dots to the position of 3C345 and 
a nearby AGN. A marginal detection was also obtained for pe- 
riod 3034 at a level of 2.1cr. It should, however, be noted that 
the small photon excess above 500 MeV has been attributed to 
a flare from Mrk 501 by Kataoka et al. (1999) because the cen- 
troid was closer to the famous TeV source, so the association of 
EGR J 1 642+3940 with 3C345 is not clear. GLAST should easily 
confirm or infirm the association. 

Several radiogalaxies (Cen A, NGC 6251, J1737-15) 
and a Seyfert 1 (GRS 1734-292) had been proposed as 
possible counterparts to 3EG sources dHartman et al., 1999, 
|Combi et al., 2003||Foschini et al., 2005l|Di Cocco et al., 2004| >' 
They triggered some interest because their identification would 
raise important questions about the origin of the y rays at large 
angle from the strongly beamed emission from the jet. We do 
not, however, confirm the spatial coincidence with EGR sources 
in the present work. All these galaxies lie well beyond the 95% 
confidence region of EGR sources. 



9. Conclusions 

We have searched for point-like sources in the reprocessed 
EGRET data from cycle 1 to 9 using new interstellar background 
models based on the most recent HI, CO, and dark gas data, as 
well as two different assumptions for the cosmic-ray distribution 
(the GALPROP diffusion model or a radial emissivity gradient 
fitted to the diffuse EGRET data). We have used the 3EG tools, 
likelihood method, procedure and significance threshold to de- 
tect sources, but have expanded the search to 3 different energy 
bands (above 100 MeV, 0.3-1 GeV, and above 1 GeV). The re- 
sulting number of detected sources has decreased by more than a 
third. Many unidentified sources, in particular among those spa- 
tially associated with the Gould Belt, are not confirmed as signif- 
icant excesses. Their emission can be explained by the additional 
interstellar emission and its structure. Several interesting coun- 
terparts to 3EG sources, such as radiogalaxies, massive binaries, 
and microquasars, are now found outside the 95% confidence 
region. We have cross-correlated the new source positions with 
large pulsar, supernova remant, pulsar wind nebulae, OB asso- 
ciations, and radiosource catalogues, yet half the sample has no 
attractive counterpart among the potential y-ray emitters. 30 new 
possible y-ray sources have also been found. 

This EGR catalog will be available in 
fits format at the Strasbourg astronomical 
Data Center (CDS) and in ASCII format at 
http://www.aim. univ-paris7.fr/EGRETj2atalogue/home. html 

Acknowledgements. We are deeply grateful to Bob Hartman for his helpful ex- 
planations about the construction of the 3EG catalogue, and to Seth Digel and 
Andy Strong for their help with the gas and Galprop maps. 
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Table A.l. The EGR catalogue 



Name RA Dec 1 

EGR J0008+7308 2.01 73.14 119.75 



b 

10.54 



EGR J0057-7839 14.46 -78.65 -57.47 



EGRJ0100+4927 15.01 49.45 124.37 



EGRJ0159-3609 29.78 -36.16 -110.71 



EGRJ0204+1505 31.00 15.09 147.75 



0.20 



EGR J0028+0457 7.06 4.95 112.15 -57.44 0.51 



EGR J0039-0945 9.75 -9.75 112.76 -72.38 0.27 



-38.47 0.53 



-13.40 0.27 



EGR J01 17+0254 19.41 2.91 135.83 -59.30 1.15 



EGRJ0141 + 1719 25.47 17.32 139.75 -43.90 0.85 



-73.04 0.59 



-44.26 0.48 



EGRJ0210-5058 32.58 -50.97 -83.84 -61.86 0.14 
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Table B.l. The EGR confused sources catalogue 



Name RA Dec 1 b 8 95 

EGRc J0225+6240 36.38 62.68 133.49 1.75 0.34 



EGRcJ08 18-46 13 124.74 -46.23 -97.25 -5.73 0.31 



EGRc J0842-4501 130.66 -45.03 -95.78 -1.68 0.26 



EGRc J0912+7146 138.15 71.77 141.44 36.44 0.62 



EGRc J0927+6054 141.91 60.91 153.55 42.15 0.67 



EGRc J1038-5724 159.61 -57.41 -74.25 0.96 0.40 



EGRc J1233-0318 188.46 -3.30 -65.69 59.28 1.04 



EGRc J1255-0404 193.78 -4.08 -55.30 58.78 0.71 



EGRc J1332-1217 203.04 -12.29 -41.69 49.36 0.56 



EGRcJ1740-2851 265.05 -28.85 -0.55 1.05 0.16 
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